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We  discuss  the  salient  features  of  a computer  program,  RangeCAD,  and  then 
translate  the  program’s  output  into  a catalogue  of  radar  cross-section  (RCS)  un- 
certainties. This  specific  catalogue  was  developed  by  NIST  to  standardize  RCS 
uncertainty  computations  at  the  various  RCS  measurement  sites.  We  check  un- 
certainty estimates  generated  by  RangeCAD  against  alternative  formulations  that 
approximate  equivalent  uncertainty  specifications.  Based  on  this  comparison,  we 
conclude  that  the  uncertainty  estimates  generated  by  RangeCAD  provide  realistic 
values  for  the  NIST  RCS  uncertainty  analysis. 

Key  words:  error  catalogue;  measurement  error;  radar;  radar  cross  section;  RCS; 
radar  signal  analysis;  uncertainty  analysis 


1.  Introduction 

The  National  Institute  of  Standards  and  Technology  (NIST)  has  been  assigned  the  task 
of  evaluating  radar  cross  section  (RCS)  measurement  uncertainty  at  various  Department 
of  Defense  sites.  As  a part  of  this  effort,  NIST  has  been  studying  a computer  program, 
RangeCAD,  that  was  developed  for  the  United  States  Air  Force  to  assist  in  making  uncer- 
tainty predictions. 

The  computer  program  RangeCAD  is  essentially  a sophisticated  worksheet  for  con- 
verting estimated  measurement-system  parameter  uncertainties  into  an  overall  estimate 
of  the  uncertainty  in  predicted  RCS  level,  using  a comprehensive  model  of  the  RCS  mea- 
surement system  for  making  these  predictions.  The  user  independently  specifies  values  for 
the  measurement  system  parameters  along  with  their  associated  variances.  Fortunately, 
as  a part  of  RangeCAD’s  development,  a fairly  extensive  error-component  database  was 
compiled,  consisting  of  measured  system  parameters  and  the  associated  variance  in  these 
measured  values.  The  existing  database  characterizes  measurements  taken  on  radars  at  the 
46th  Test  Group’s  RATSCAT  division,  Holloman  Air  Force  Base,  and  at  the  contractor’s 
Electronic  Test  Site.  A sample  plan  for  taking  measurements  to  update  this  database  or 
to  customize  it  for  other  RCS  facilities  is  described  in  the  program  users  manual. 

RangeCAD  is  capable  of  fully  accounting  for  all  target-ground  reflections  occur- 
ring within  the  specified  range-gate  window,  incorporating  the  effects  of  uncertainties 
attributable  to  the  transmitted  waveform,  the  antenna  pattern,  and  the  receiver  in  its 
predictions.  It  can  also  provide  a premeasurement  diagnosis  of  the  measurement  setup. 
The  program’s  principal  limitation  is  that  a model  of  the  target’s  radar  cross  section  is 
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needed  in  order  to  determine  the  measured  RCS  of  the  target.  This  is  handled  by  utilizing 
the  program  to  construct  a dynamic-range  model  of  the  measurement  system,  from  which 
the  expected  mean  and  variance  of  the  measured  RCS  can  be  predicted  for  any  measure- 
ment level.  Construction  of  the  dynamic-range  model  is  accomplished  by  replacing  the 
target  with  a single  point  scatterer  whose  scattering  amplitude  varies  over  the  expected 
measurement  range.  The  field  taper  across  the  spatial  extent  of  the  target  remains  un- 
accounted for  in  this  implementation.  However,  the  program  permits  field-taper  mapping 
by  displacing  the  point  scatterer  from  the  center  of  the  target  mount,  thereby  enabling 
a dynamic-range  model  to  be  produced  for  any  desired  location  within  the  target  zone’s 
spatial  extent. 

The  output  from  RangeCAD  can  be  incorporated  into  the  NIST  RCS  uncertainty 
analysis,  which  is  a table  of  significant  error  sources  designed  to  facilitate  insight  into 
RCS  uncertainty  computations  and  to  standardize  these  calculations.  The  estimates  of 
RCS  uncertainty  generated  by  RangeCAD  can  be  compared  to  uncertainty  estimates  ob- 
tained from  alternative  formulations  (see  below).  In  spite  of  differences  between  underlying 
assumptions,  the  results  generated  by  RangeCAD  and  the  alternative  formulations  were 
consistent  within  an  acceptable  margin,  an  attribute  which  tended  to  lend  credibility  to 
RangeCAD. 


2.  Salient  Features  of  RangeCAD 


RangeCAD  is  designed  to  model  a ground-bounce  RCS  antenna  range  for  targets  consist- 
ing of  an  ensemble  of  spatially-distributed  scattering  centers  (point  scatterers).  Use  of 
such  models  to  simulate  extended  targets  has  proven  very  successful  [1].  The  RangeCAD 
program  fully  accounts  for  direct  and  ground-reflection  ray  paths  from  the  antenna  to  each 
point  scatterer  on  the  target  and  back  again.  Also,  in  the  case  of  targets  specified  using 
GTD  (Geometrical  Theory  of  Diffraction),  coupling  between  the  target  and  the  ground 
is  accounted  for  by  adjusting  the  target’s  scattering-point  locations  to  correspond  to  a 
ground-bounce  ray-path  excitation.  The  round-trip  propagation-path  length  is  obtained 
for  each  scattering  point  along  all  ray  paths  in  order  to  determine  the  corresponding  exci- 
tation amplitudes  on  the  transmitted  waveform’s  envelope.  The  program  determines  the 
received  voltage  from  each  scattering  point  as  a function  of  the  transmitted  waveform  am- 
plitude, the  antenna  gain  in  the  direction  of  the  scatterer,  the  range  propagation  gain,  the 
reflected  amplitude  and  phase  of  the  scatterer,  and  coherent-integration  Doppler-mismatch 
loss.  The  instantaneous  received  voltages  from  each  point  scatterer  within  the  receiver’s 
range-gate  window  are  summed  along  with  the  background  clutter,  and  then  a receiver 
linearity  correction  is  applied  to  obtain  the  detected  signal.  A separate  calculation  of 
additive  noise  is  made  to  include  receiver  noise,  transmitter  noise,  antenna  noise,  and  rf 
interference. 

The  instantaneous  received  power  from  a single  scatterer  [2]  is  given  by 


PR{t)  = PT(t  - r) 


^ ^ 

(47r) 


(1) 


where  r = 2r/c  is  the  scatterer’s  round  trip  delay  time,  r is  the  range  to  the  scatterer,  c 
is  the  propagation  velocity,  A is  the  wavelength,  G is  the  one-way  antenna  power  gain  in 
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the  direction  of  the  scatterer,  Pt(0  instantaneous  transmitted  power,  and  a is  the 

scatterer’s  RCS.  A ground-reflection  factor  must  be  inserted  when  the  propagation  path 
intersects  the  ground.  The  instantaneous  receiver  response  (video  signal)  to  the  point 
scatterer  is  given  [2]  (also  see  Appendix  A)  by 


= x(^  - 'r)D{u) 


G\ 

(47r)tr^ 


i 
(7  2 


(2) 


which  has  a structure  similar  to  eq  (1).  Here,  D{u)  represents  Doppler  integration- 
mismatch  loss  for  a Doppler-shift  frequency  v,  while  x(0  ambiguity  function)  repre- 
sents a convolution  between  the  receiver  impulse  response  and  the  pulse  shape  (envelope) 
of  the  transmitted  signal.  Such  a convolution  is  included  in  RangeCAD’s  waveform  model. 
For  targets  consisting  of  an  ensemble  of  scatterers,  the  individual  complex  video  signals 
are  summed  [2].  After  a receiver  linearity  correction  is  made  in  RangeCAD,  the  resulting 
sum  is  squared  and  then  combined  with  an  additive  noise  term  to  obtain  the  composite 
video  power  sampled  by  the  receiver  [2,3]. 

In  RangeCAD,  the  power  received  from  a target  is  calculated  as  a function  of  that 
target’s  ideal  RCS.  The  measured  RCS  of  a target  is  obtained  from  the  ratio  of  the  target’s 
received  power  to  the  calibration  standard’s  received  power  times  the  calibration  standard’s 
ideal  (theoretical)  RCS  cri^\  The  measurement  error  is  then  obtained  as  the  difference 
between  the  target’s  ideal  and  measured  RCS.  The  measured  RCS  (Tm  is  given  by 


<^m  — 


(3) 


Here,  Pt  and  Pg  represent  calculated  received  power  from  the  target  and  the  standard, 
respectively,  based  on  a simulation  of  the  measiorement  system  (including  system  errors). 

Program  RangeCAD  supports  ideal  RCS  calculations  for  three  built-in  target  models. 
For  two  of  these  models,  a thin  flat  plate  and  a finite-length  circular  cylinder,  the  ideal 
RCS  is  calculated  using  GTD  analysis.  For  the  third  model,  a sphere,  the  ideal  RCS  is 
calculated  using  the  Mie  series.  In  addition  to  these  theoretical  models,  the  ideal  RCS  for 
an  arbitrary  target  may  be  specified  by  an  array  of  point  scatterers,  in  which  case  the  user 
must  specify  the  scattering  amplitude,  phase,  and  position  of  each  point  scatterer. 

A single  point  scatterer  may  be  specified  to  simulate  a probe  target.  For  more  complex 
targets,  a subsidiary  program  is  available  to  generate  a planar  array  of  point  scatterers 
corresponding  to  a given  far-field  pattern,  such  as  a computed  scattered  field.  In  principle, 
a planar  array  of  point  scatterers  could  be  generated  from  measured  bistatic  RCS  data 
to  represent  an  arbitrary  target.  RangeCAD  would  then  be  able  to  give  a first-order 
approximation  to  the  measurement  uncertainty.  Generally,  such  a procedure  is  impractical 
due  to  the  excessive  amount  of  data  required. 

As  a consequence,  the  normal  use  of  RangeCAD  is  to  obtain  a family  of  simulated 
RCS  measurements  from  a single  point-scatterer  target  whose  scattering  amplitude  varies 
over  the  measurement  system’s  dynamic  range.  Thus,  RangeCAD’s  primary  usefulness 
lies  in  enabling  the  user  to  evaluate  the  uncertainty  in  an  actual  RCS  measurement  from 
a lookup  table  incorporating  a dynamic-range  model  of  the  measurement  system.  In  this 
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role,  RangeCAD  can  evaluate  the  effect  of  system  parameter  changes  on  RCS  measurement 
systems  and  provide  premeasurement  diagnoses.  A separate  postmeasurement  processing 
program  should  be  used  to  compensate  actual  RCS  measurements  for  recognized  systematic 
errors,  such  as  receiver  nonlinearity,  and  the  benefit  of  these  corrections  incorporated  into 
RangeCAD’s  input  parameter  set. 

The  uncertainty  in  a RCS  measurement  can  be  expressed  in  decibels  as 


101og(l  ± 


(4) 


where  /ig.  is  the  total  RCS  variance  and  a is  the  cissociated  RCS  level  (nominally,  a either 
denotes  the  measured  RCS  of  the  calibration  standard  or  the  target).  If  the  RCS  is  a 
function  of  a set  of  N parameters  Xi,  and  fix,  is  the  variance  of  r,-,  i = 1,  iV,  then  the  total 
RCS  variance  [4]  is  given  by 


N 


= E! 

i=l 


da 

dxi 


(5) 


RangeCAD  only  includes  those  parameters  in  the  above  sum  which  have  been  assigned 
nontrivial  variance  values  by  the  user.  In  eq  (5),  the  coefficient  of  Hxi  is  called  the  sensi- 
tivity, and  the  product  of  that  sensitivity  with  fixi  represents  a partial  RCS  variance  due 
to  the  uncertainty  in  Xi  . The  partial  derivatives  are  calculated  using  a simple  two-term, 
divided- difference  formula.  That  is,  the  derivative  is  obtained  by  calculating  the  RCS  at 
a small  increment  above  and  below  the  mean  value  for  the  parameter,  and  dividing  the 
difference  between  these  two  results  by  two  times  that  increment.  In  terms  of  Lagrange 
interpolation,  this  approximation  to  the  first  derivative  is  correct  through  to  the  third 
derivative  [5].  As  a result,  whatever  RCS  functional  dependence  a user  might  specify  for 
a parameter,  that  dependence  is  fully  taken  into  account  by  the  sensitivity  calculation. 

RangeCAD  was  designed  so  that  separate  data-input  files  are  specified  for  the  measure- 
ment of  calibration  standards  and  the  measurement  of  targets.  This  feature  is  normally 
used  to  change  the  rf  attenuation  between  the  cahbration  measurement  and  the  target 
measurement.  Radar  parameter  measurements  have  been  carried  out  at  VHF,  UHF,  L,  S, 
C,  and  X band  frequencies,  and  corresponding  RangeCAD  data  files  have  been  developed 
[6]  for  use  at  the  contractor’s  Electronic  Test  Site  and  at  RATSCAT.  Generally,  mean 
parameter  values  are  given,  while  the  variances  represent  the  measurement  uncertainty. 
Only  the  variances  of  those  parameters  which  resulted  in  significant  RCS  uncertainties  are 
retained  in  the  data  files. 

Although  RangeCAD  was  originally  designed  for  use  at  RATSCAT,  the  program’s 
modular  nature  allows  it  to  be  readily  adapted  for  use  at  other  sites.  For  instance,  the 
U.S.  Navy’s  Santa  Cruz  Radar  Imaging  Facility  (SCRIF)  has  the  same  model  radar  ais  at 
RATSCAT,  so  RangeCAD  could  be  put  to  use  at  SCRIF  using  very  similar  input  data 
as  that  generated  at  RATSCAT.  Step-by-step  procedmres  for  carrying  out  and  verifying 
parameter  measurements  are  given  in  the  RangeCAD  user’s  manual  [6],  so  the  database 
can  readily  be  customized.  For  instance,  the  basic  parameter  set  has  been  modified  at 
RATSCAT  to  characterize  three  different  measurement  systems.  A significant  feature 
enabling  RangeCAD  to  be  used  at  other  sites  is  the  separate  data-input  file  structure 
mentioned  above.  For  instance,  at  SCRIF  the  calibration  standard  is  a dihedral  located 
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across  a valley,  whereas  the  targets  are  ships,  so  for  this  facility  the  calibration  measure- 
ment’s ground-plane  parameters  woiild  be  set  to  free-space  values,  while  for  the  target’s 
measurement  they  would  simulate  the  sea.^ 

3.  RangeCAD  and  the  NIST  RCS  Uncertainty  Analysis 

The  NIST  uncertainty  analysis  for  RCS  measurements  was  created  to  identify  imcertainty 
categories  that  would  apply  throughout  the  RCS  measurement  community  [7].  Subsidiary 
uncertainty  tables  specific  to  individual  sites  can  easily  be  incorporated  within  this  frame- 
work. Table  1 presents  an  example  of  the  NIST  RCS  uncertainty  analysis.  Here,  separate 
tables  are  provided  for  the  target  and  calibration  measurements.  Partial  RCS  uncertain- 
ties are  combined  on  the  last  line  of  each  table  by  converting  the  partial  uncertainties 
(expressed  in  decibels)  to  variances  and  then  summing  these  variances  to  obtain  the  total 
RCS  variance.  The  overall  RCS  uncertainty  in  decibels  is  then  obtained  from  the  total 
RCS  variance  as  shown  in  eq  (4). 

Predicted  RCS  uncertainties,  obtained  from  RangeCAD  using  an  S-band  data  set, 
are  incorporated  into  the  NIST  RCS  uncertainty  analysis  framework  as  shown  in  Table 
1.  The  numerical  values  presented  in  the  table  should  be  considered  illustrative  rather 
than  applicable  to  any  specific  measurement.  Moreover,  the  computed  uncertainty  values 
depend  upon  the  RCS  level.  In  the  case  of  Table  1,  the  values  presented  assume  distinct 
calibration  and  target  RCS  levels  of  10  dBsm  and  —20  dBsm,  respectively.  Obviously,  a 
number  of  tables  would  be  needed  to  fully  characterize  a RCS  measurement.  The  uncer- 
tainty associated  with  specific  RCS  measurement  levels  could  then  be  determined  from 
the  appropriate  table.  Alternately,  a subsidiary  program  is  available  with  RangeCAD  for 
automating  this  process. 

In  the  following,  each  entry  in  the  NIST  RCS  uncertainty-analysis  table  is  examined, 
and  RangeCAD ’s  capability  for  estimating  the  partial  RCS  uncertainty  is  indicated  where 
applicable.  Alternative  formulations  for  validating  specific  results  axe  also  discussed. 

1.1  & 2.1  Average  Illumination:  Tracking  and  pointing  uncertainty  constitutes 

this  error  category  [7].  Specifying  allowable  variances  in  the  target’s  position  variables 
enables  RangeCAD  to  calculate  the  uncertainty  directly.  Alternately,  if  an  angular  tracking 
uncertainty  A is  specified,  then  (see  Appendix  B)  the  RCS  uncertainty  in  decibels  is 
— 201og[l  — sin^(  jA/^o)  sinc(2X)].  Here,  sincX  = sinX/X,  where  X = is  the 

angular  span  of  the  target,  and  is  the  half-power  beamwidth.  Alternately,  we  can  obtain 
X from  a measurement  of  the  FieldTaper  (cf.  Near  Field  below).  Our  calculations  in  the 
table  use  the  values  FieldT aper  = 0.5  dB,  = 2.5°,  and  A = 0.1°. 

1.2  & 2.2  Background/Target  Interactions:  Target/pylon  interaction  and  clutter 
data  are  supplied  to  RangeCAD  by  files  TGT_MOUNT.DAT  and  CAL_MOUNT.DAT, 
which  contain  tabular  frequency  and/or  azimuth  dependent  values  which  can  be  used  for 
interpolation.  Target-to-pylon  RCS  ratios  as  small  as  10  dB  can  result  in  a measurement 

^ The  prospective  user  should  verifythat  a program  modification  suggested  by  NIST  has 
been  implemented  so  that  the  free-space  ground-plane  response  will  be  correct. 
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uncertainty  of  3 dB.  As  the  target  RCS  approaches  the  pylon  RCS  destructive  interference 
effects  can  be  dramatic  [8].  Numerical  data  for  specifying  file  values  are  not  available  to 
NIST. 

1.3  & 2.3  Cross-Polarization  Effects:  Measurements  carried  out  [9]  at  RATSCAT 
in  1988  indicated  a cross- polarization  isolation  of  less  than  -30  dB.  The  values  in  the  table 
assume  a cross-polarization  level  around  -40  dB  with  a non-depolarizing  standard  and  a 
depolarizing  target.  RangeCAD  may  be  used  directly  to  obtain  the  uncertainty  by  spec- 
ifying a polarization- angle  variance  in  RangeCAD ’s  antenna-parameter  file.  Alternately 
[7],  as  in  the  table,  assuming  a cross-polarization  level  of  — e dB  and  a depolarizing  target 
produces  an  uncertainty  of  201og(l  -1-2  x 10“^/^°),  while  a non- depolarizing  target  produces 
an  uncertainty  of  201og(l  -f  10“^/^®). 

1.4  & 2.4  Drift:  Tests  at  RATSCAT  carried  out  in  1991  [6]  and  in  1988  [9]  indicated 
propagation  path  drifts  of  about  0.5  dB  over  an  hour’s  time.  Periodic  recalibration  us- 
ing secondary  standards  limits  expected  drift.  The  table  assumes  neghgible  drift  during 
calibration. 


1.5  &:  2.5  Frequency:  As  discussed  above,  RangeCAD  obtains  the  RCS  sensitiv- 
ity to  a parameter  by  recalculating  the  RCS  at  small  increments  above  and  below  the 
parameter’s  mean  value.  For  RCS  frequency  response,  RangeCAD  supports  specifying 
target-pylon  clutter  as  a function  of  frequency,  thereby  explicitly  enabling  target/target- 
mount  frequency  effects  to  be  modeled.  Inasmuch  cls  frequency  variance  is  Hmited  to  the 
bandwidth  squared,  except  for  unusual  causes  the  expected  contribution  to  RCS  variance 
from  gradients  in  the  received- signal  amplitude  as  a function  of  frequency  should  be  small. 
However,  a small  change  in  frequency  can  result  in  a significant  phase  change  between  scat- 
terers,  which  could  lead  to  a large  RCS  variance.  For  example,  two  scatterers  separated 
in  range  by  L will  undergo  a 27r  relative  phase  change  in  RCS  due  to  a frequency  change 
A/,  where  A///  = X/{2L)  [2].  For  T/A  = 125,  a 0.2%  relative  frequency  shift  results  in 
a 180°  phase  shift,  which  could  cause  the  combined  response  to  vary  from  a maximum  to 
a null. 


1.6  2.6  Integration:  Both  A/D  quantization  loss  and  Doppler  integration  mismatch 

loss  fit  this  category.  Integration  uncertainty  is  obtained  by  summing  the  individual  par- 
tial RCS  variances  and  using  eq  (4)  to  convert  the  result  to  a partial  RCS  uncertainty 
expressed  in  decibels.  Assuming  an  8-bit  A/D  converter,  the  quantization  loss  is  given  by 
— 20  log(255/256)  0.03  dB.  Doppler  integration-mismatch  loss  is  given  [6]  (cf.  (A-11)) 
by 


Nsmy 


The  pulse  phase  shift  (j)  = 27r  x width  x AzimuthRotationRate / {\  x PRF),  where  width 
denotes  the  target’s  (electrical)  span,  A the  wavelength,  PRF  the  pulse  repetition  fre- 
quency, and  N the  number  of  integrated  pulses.  Normally,  Doppler  integration-mismatch 
loss  can  be  expected  to  be  small. 


1.7  h 2.7  I/Q  Imbalance:  The  table  directly  uses  RangeCAD’s  predictions.  Target 
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and  calibration-standard  uncertainty  calculations  are  independently  carried  out  at  distinct 
RCS  levels.  Gain  imbalance  and  phase  imbalance  contribute  to  this  category.  The  indi- 
vidual partial  RCS  variances  are  summed  and  then  eq  (4)  is  used  to  convert  the  result  to 
a partial  RCS  uncertainty  expressed  in  decibels. 

1.8  & 2.8  Near  Field:  Field  taper  across  the  target  constitutes  this  error  category. 

Uncertainty  may  be  obtained  directly  from  RangeCAD  using  the  program’s  probe  option 
to  predict  RCS  gain  reduction  in  the  target  volume.  Alternately,  uncertainty  can  be  de- 
termined from  geometry  and  the  antenna’s  half-power  beamwidth  9q  or  from  field-taper 
measurements.  For  instance,  probe  measurements  [9]  in  1988  showed  horizontal  and  ver- 
tical field  tapers  of  0.5  dB.  Transmitting  antenna  gain  and  receiving  antenna  gain  doubles 
the  potential  error.  Assuming  a cosine- squared  antenna  pattern  and  a target  that  is  a lin- 
ear array  of  equal-amphtude  uniformly-distributed  point  scatterers,  the  RCS  imcertainty 
[10]  is  — 201og[(l-f-sinc(2X))/2],  where  sincJA  = sinX/X.  Here,  X = where  9i 

is  the  angle  subtended  by  the  target  at  the  antenna.  Alternately,  one  may  obtain  X by 
solving  the  equation  — 201og(cos  AT)  = FieldTaper.  The  latter  option  is  used  in  the  table, 
where  FieldTaper  = 0.5  dB,  in  accord  with  measured  values  reported  above. 

1.9  &:  2.9  Noise/ Background:  The  table  directly  uses  RangeCAD’s  predictions. 
Target  and  calibration-standard  uncertainty  calculations  are  independently  carried  out  at 
distinct  RCS  levels.  Noise,  background  magnitude,  clutter  amplitude,  backgroimd  phase, 
clutter  phase,  and  ground-plane  roughness  all  contribute  to  this  category.  The  partial  RCS 
variances  are  summed  and  then  eq  (4)  is  used  to  obtain  this  category’s  combined  partial 
RCS  uncertainty. 

1.10  Sz  2.10  Nonlinearity:  RangeCAD  predictions  are  used  directly  in  the  table. 
Target  and  calibration-standard  uncertainty  calculations  are  independently  carried  out  at 
distinct  RCS  levels.  Receiver  nonlinearity  variance  and  a number  of  waveform  parameter 
variances  all  contribute  to  this  category.  The  individual  partial  RCS  variances  are  siunmed 
and  then  eq  (4)  is  used  to  obtain  the  combined  partial  RCS  imcertainty  for  this  category. 
Inasmuch  as  a point-scatterer  target  is  used,  the  waveform-parameter  variances  should  be 
set  equal  to  the  square  of  the  maximum  variation  that  a point  scatterer  would  experience 
on  a range  walk  over  the  extent  of  the  actual  target. 

1.11  & 2.11  Range:  Assuming  that  the  target  range  uncertainty  is  AR,  the  RCS 
uncertainty  is  — 401og(l  — AR/R),  which  is  usually  negligible.  RangeCAD  may  be  used 
directly  for  this  calculation. 

1.12  & 2.12  Target  Orientation:  Uncertainty  may  be  obtained  from  RangeCAD  by 
specifying  allowable  variances  in  the  target’s  angular  positioning  variables.  Normally,  this 
uncertainty  contribution  can  be  expected  to  be  small. 

2.13  Computational  Error:  Numerical  computation  error  for  standard  target’s 
RCS,  including  manufacturing  defects.  The  table’s  estimate  of  0.5  dB  is  consistent  with 
RCS  uncertainty  estimates  [9]  for  cylinders  and  with  uncertainties  in  numerical  computa- 
tions obtained  at  NIST. 

1.13  Calibration  Uncertainty:  Overall  calibration  uncertainty  from  2.14. 
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1.14  Sz  2.14  Overall  Uncertainty:  The  partial  RCS  variances  are  obtained  for  each 
category  and  then  summed  to  obtain  the  total  RCS  variance.  The  overall  imcertainty 
expressed  in  decibels  is  then  obtained  using  eq  (4). 

4.  Conclusion 

Program  RangeCAD’s  formulation  for  the  radar  signal  received  from  an  illuminated  target 
agrees  in  form  with  established  expressions  for  representing  the  received  video  signal  of  a 
monostatic  radar  due  to  rf  illumination  of  an  ensemble  of  point  scatterers.  We  have  pointed 
out  significant  features  of  RangeCAD  and  discussed  the  possibihty  of  using  program  out- 
put to  ease  the  experimentalist’s  task  of  constructing  an  RCS  uncertainty  analysis.  In 
this  regard,  when  using  dynamic-range  modeling  to  evaluate  RCS  uncertainty  it  is  appar- 
ent that  some  partial  variances  should  be  exaggerated  (for  example,  waveform  parameter 
variances  should  be  augmented  to  include  the  range  of  waveform  values  experienced  over 
the  downrange  extent  of  the  target).  Also,  some  variances  (for  example,  near  field  or  field 
taper)  are  readily  evaluated  by  alternative  expressions  whose  use  may  be  more  convenient. 
The  RangeCAD  program  is  applicable  to  or  can  be  adapted  to  a wide  variety  of  RCS 
ranges.  Although  RangeCAD  is  not  a panacea  for  determining  RCS  uncertainty,  it  can  be 
viewed  as  a useful  tool  to  steer  the  experimentalist  into  considerations  that  might  otherwise 
be  overlooked,  and  it  could  become  a model  for  producing  comparable  RCS  uncertainty 
estimates  at  applicable  ranges.  Moreover,  the  task  of  constructing  the  RangeCAD  pa- 
rameter tables  enables  the  user  to  track  equipment  stability  over  time,  and  could  lead  to 
RCS  measurement  error  detection  (as  opposed  to  RCS  imcertainties)  which  would  require 
compensation  by  a post-measurement  processing  program.  Finally,  using  RangeCAD  to 
carry  out  premeasurement  diagnosis  might  avoid  wasting  valuable  measurement  time  by 
detecting  otherwise  overlooked  problems. 


The  authors  would  like  to  express  their  appreciation  to  Barry  Stroman,  Martin  Marietta 
Orlando  Aerospace,  for  helpful  discussions  regarding  RangeCAD. 
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Table  1.  Typical  RCS  uncertainties  based  in  part  on  RangeCAD  data 


MEASUREMENT  UNCERTAINTIES 

+dB 

-dB 

1.1 

Average  illumination 

0.03 

-0.03 

1.2 

Background /target  interaction 

? 

? 

1.3 

Cross  polarization 

0.16 

-0.17 

1.4 

Drift 

0.50 

-0.57 

1.5 

Frequency 

0.00 

0.00 

1.6 

Integration 

0.03 

-0.03 

1.7 

I/Q  imbalance 

0.67 

-0.79 

1.8 

Near  field 

0.33 

-0.35 

1.9 

Noise/background 

0.65 

-0.76 

1.10 

Nonlinearity 

0.07 

-0.07 

1.11 

Range 

0.00 

0.00 

1.12 

Target  orientation 

0.00 

0.00 

1.13 

Calibration  uncertainty  (2.14) 

0.83 

-1.02 

1.14 

Overall  uncertainty 

1.29 

-1.85 

CALIBRATION  UNCERTAINTIES 

+dB 

-dB 

1.1 

Average  Illumination 

0.03 

-0.03 

1.2 

Background /target  interaction 

? 

? 

1.3 

Cross  polarization 

0.00 

0.00 

1.4 

Drift 

0.00 

0.00 

1.5 

Frequency 

0.00 

0.00 

1.6 

Integration 

0.03 

-0.03 

1.7 

I/Q  imbalance 

0.66 

-0.77 

1.8 

Near  field 

0.16 

-0.17 

1.9 

N oise /background 

0.12 

-0.13 

1.10 

Nonlinearity 

0.05 

-0.05 

1.11 

Range 

0.00 

0.00 

1.12 

Target  orientation 

0.00 

0.00 

1.13 

Computational  error 

0.50 

-0.57 

1.14 

Overall  calibration  uncertainty 

0.83 

-1.02 

10 


Appendix  A.  Radar  Signal  Representations 


Here  we  present  a brief  review  of  the  mathematical  representation  of  the  signal  received  by 
a radar  and  the  associated  receiver  response.  Additional  details  may  be  found  in  [2].  The 
theory’s  application  to  waveform  analysis  and  range-Doppler  resolution  is  further  discussed 
in  [3]. 

Let  tprii)  denote  a transmitted  rf  signal  at  time  t having  power  Pt(0  — 
let  7 denote  the  complex  reflection  coeflacient  of  a point  target  (i.e.,  a target  that  is  small 
compared  to  a resolution  cell)  such  that  the  point  target’s  radar  cross  section  a = |7p. 
Then  the  signal  received  by  the  radar  from  the  scatterer  7 is 


- r{t)) 


GX 


(4'7r)  2 r' 


:7, 


(.4-1) 


where  A is  the  wavelength,  r is  the  range,  G is  the  one-way  antenna  gain  in  the  direction 
of  the  scatterer,  and  r(t)  is  the  round-trip  delay.  The  received  power  corresponding  to  the 
received  signed  ipR{t)  results  in  the  familiar  radar  range  equation 


PR{t)  = \^R{t)\'^  = Prit  - r(t)) 


(47r)^r^ 


(A -2) 


We  now  consider  that  the  complex  transmitted  signal  can  be  decomposed  into 

a sinusoided  function  of  the  carrier  frequency  / = c/A,  where  c is  the  speed  of  light,  and 
a complex  modulation  function  containing  only  low  frequencies,  as 

We  further  consider  that  the  target  of  interest  consists  of  an  ensemble  of  point  scatterers. 
Modeling  extended  targets  as  a collection  of  point  scattering  centers  is  discussed  in  [1].  The 
received  signal  from  the  target  is  then  given  as  a superposition  of  the  component  signals 
from  each  scattering  center.  We  let  jk,  Gkt  and  rk  denote  the  scatterer’s  complex 
reflection  coefficient,  gain,  and  range,  respectively.  The  instantaneous  round-trip  delay 
can  be  expanded  into  a Taylor  series  as  r(t)  = rjt  — Ukt/f^  where  Tk  = 2rk/c  is  the  delay 
for  the  scatterer  and  Uk  = — 2f/A  is  the  associated  Doppler  coefficient  (this  assumes 
the  scatterer  to  be  slowly  traveling  at  a constant  range  rate  r over  some  time  interval  of 
interest).  Combining  these  specifications  with  our  previous  expression  (A-1)  for  the  signal 
received  from  a single  point  scatterer,  we  obtain  [2] 


Mt)  = 


(A -3) 


where 


A 


(47r)t 


(A -4) 


We  now  consider  that  the  transmitted  signal’s  complex  modulation  function  consists 
of  a train  of  M short  pulses,  so  that 


M-l 


^ - nTr), 


(A -5) 


n=0 
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where  Tr  is  the  pulse  repetition  period.  Now  if  the  individual  pulse  length  Tp  is  less  than 
half  Tr,  |rjfc|  <C  MTr,  and  \i/k\  then  we  can  introduce  a much  simplified  analysis 

(thereby  allowing  us  to  evaluate  each  pulse’s  response  cls  equivalent  to  that  from  a central 
pulse).  Let  rj.  = rfc(mod  Tr)  and  = i/fc(mod  l/T^).  Consequently,  0 < rj.  < Tr  and 
0 < \u^\  < 1/Tr.  Now  for  the  moment,  let  us  define  t as  the  time  measured  from  the 
transmission  of  the  pulse  and  confine  t as  0 < t < Tr.  Then  we  can  replace  t in  our 
expression  ( A-3)  for  the  received  signal  by  f + nTr  and  write  the  received  signal  due  to  the 
pulse  as 

+ nTr)  = ^ Vkfippit  - {A  - 6) 

it 

where  we  have  dropped  the  essentially  constant  phase  term  2TTi''f.t  (having  bounded  vari- 
ation 27r  1 I Tp  <C  27r).  The  received  signal  is  usually  mixed  down  to  video  as  the  initial 
processing  step  in  the  receiver,  so  the  received  video  signal  can  be  written  as 

^iR{tAntr)  = '^Vkfiprit  - {A  - 7) 

it 

If  pulse  compression  is  implemented  next,  then  the  output  [2]  corresponding  to  the 
pulse  is  given  by 

T 

ZA-r)  = / " mit  + nTr},x-pF{t  - T)dt  = V VkXp{r  - (A  - 8) 

Jo  i 

where  fJ^ppi—t)  is  the  impulse  response  of  the  receiver’s  filter  (the  receiver  is  assumed  to  be 
matched  to  a pulse-signal  ^ippit)  that  resembles  the  pulse  fxpxii)  that  is  transmitted  by  the 
radar,  and  the  asterisk  denotes  complex  conjugation).  Under  our  short-pulse  assumption 
the  filter  response  or  ambiguity  function  XP{'^)  becomes  independent  of  Doppler  shift  and 
reduces  to  the  cross  correlation 

/-l-oo 

fiPT{t)fi*pF{-t  - r)dt.  {A  - 9) 

-OO 

Next,  if  coherent  integration  of  N pulses  is  implemented,  N < M,  then  we  would 
sum  the  individual  pulse  responses  given  in  eq  (A-8)  over  n.  We  note  that  at  this  point  a 
tuning  frequency  could  be  introduced  to  compensate  the  received  signal  for  Doppler  shift; 
however,  here  we  assume  that  the  receiver  is  kept  tuned  to  the  transmitted  frequency, 
which  is  applicable  to  airborne  target  simulation.  If  we  also  assume  square-law  detection, 
our  final  expression  for  the  detected  signal  becomes  [2,3] 


iV-l 


\Z{r)f 


= = |^Uxp(r  - (A  - 10) 


n=0 


It  may  be  further  noted 
D(u)  can  be  written  as 


that  the  absolute  value  of  the  Doppler  fine  structure  coefficient 


sm{N7cvTr) 

N sin(7ri/Tr) 


(A -11) 
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Appendix  B.  Estimation  of  Amplitude  Taper  Uncertainty 


Here  we  present  a brief  review  of  the  estimated  RCS  uncertainty  due  to  amplitude  taper 
across  a long  target.  In  order  to  make  the  analysis  tractable,  a number  of  simplifying 
assumptions  are  made  to  characterize  a generic  RCS  target.  In  particular,  the  target 
is  modeled  as  a uniform  distribution  of  equal  amplitude  point  scatterers.  A nonuniform 
distribution  function  can  be  used  when  apriori  information  is  available  about  the  target. 
Additional  details  concerning  this  analysis  may  be  found  in  ref.  [10]. 

The  positions  of  the  scatterers  along  the  transverse  dimensions  of  the  target  are  de- 
scribed by  a probability  density  function  p{6,  (f))  with  the  radar  located  at  the  origin  of  a 
standard  spherical  coordinate  system.  The  ratio  of  the  measured  RCS  to  the  standard’s 
RCS  is  given  by 


3^2 

^R{t) 


(B-1) 


Let  the  reference  target  consist  of  a single  point  scatterer  70.  Then  the  ratio  of  received 
signals  can  be  written  as 


V’r(<)  l^rit  - to)Go7o 


(B-2) 


In  most  cases,  the  amplitude  variations  caused  by  differences  in  range  between  the  scat- 
terers is  negligible.  We  further  suppose  that  we  can  combine  the  amplitude  modulation, 
the  phase  discrepancy  between  each  scatterer,  and  the  scatterer’s  reflection  coefficient  to 
produce  an  expected  real  reflection  coefficient  7^.  We  further  assume  that  the  gain  in  the 
direction  of  the  calibration  target  is  unity.  With  these  simplifications  we  can  write  the 
expected  absolute  value  ratio  of  the  received  signals  as 


^{1 


v-g(<) 

>/’»(<) 


70  V 


(R-3) 


If  we  further  assume  all  scatterer’s  have  the  same  cross-section  and  all  of  their  positions 
are  described  by  p(^,  <j))^  then  the  expected  absolute  value  ratio  becomes 

^ { I I } - / f p{e,4>)G{e,<i>)ded4,.  (b  - 4) 

7o  J J 

In  order  to  estimate  the  effects  of  amplitude  taper,  we  will  form  the  ratio  of  the  measured 
signal  to  the  signal  obtained  when  the  target  is  in  the  far  field.  We  also  assume  a long  thin 
target  with  the  long  dimension  transverse  to  antenna  boresight.  Consequently,  the  double 
integration  can  be  transformed  to  a single  integral.  Using  the  variable  y = r9,  where  r is 
the  range  to  the  target,  we  obtain 

I !"///!!  I } - / i^-^) 

Wr  (u 
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where  now  p{y)  describes  the  point-scatterer  probability  density  along  the  target  length 
L.  The  ratio  7^/70  was  eliminated  by  forming  the  ratio  of  the  measured  signal  to  the 
predicted  far-field  response  [10]. 

Now  if  the  center  of  the  target  is  offset  by  8y  from  antenna  boresight  and  we  assume 
a uniform  point-scatterer  probability  density  function,  then 


^{1 


n(t)  p ^ 1 

^ J-k+h 


G{y)dy. 


{B-6) 


A simple  antenna  gain  model  that  is  easy  to  integrate  is  the  cosine  squared  pattern.  Ac- 
cordingly, we  assume  G{rd)  = cos^(-|^),  where  is  the  antenna’s  half-power  beamwidth. 
If  we  also  assume  Sy  = rA,  where  A is  the  angular  tracking  error,  and  we  define  L = rBi, 
where  61  is  the  angle  sustended  by  the  target,  then  we  obtain 


^{i 


II  ^ 1 

’I’ltit)  i-lSi+A 


o.TT  6 . 


Carrying  out  the  integration  we  readily  obtain 


(B-7) 


V’ff  (0  n ^ r-,  • ^ ^ M 


(B-8) 


where  we  have  made  use  of  the  well  known  definition,  sine  A"  = sin  A/ A. 

The  preceding  expression  (B-8)  gives  the  effects  of  both  amplitude  taper  across  the 
target  and  pointing  error.  In  carrying  out  an  uncertainty  analysis,  we  separate  these 
two  effects.  We  set  A = 0 in  eq  (B-8),  and  express  in  decibels  the  total  RCS  near-field 
uncertainty  as 


20  log 


2 

1 +smc(f|j)' 


(B-9) 


We  can  also  express  in  decibels  the  pointing  uncertainty  as 


20  log  ( 1 + sinc( 


{B  - 10) 


This  expression  was  obtained  by  subtracting  the  expression  for  the  effects  of  both  amplitude 
taper  and  pointing  error  from  the  expression  for  amplitude  taper  alone  and  adding  unity. 
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